function [C]=ShannonCost(mu,p,gamma_g,gamma_b)

C=-1*(H(mu)+ H(1-mu)- p*(H(gamma_g)+H(1-gamma_g))-(1-p)*(H(gamma_b)+H(1-gamma_b)));
end

